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We observe chimeralike states in an ensemble of oscillators using a type of global coupling consist¬ 
ing of two components: attractive and repulsive mean-field feedback. We identify existence of two 
types of chimeralike states in a bistable Lienard system; in one type, both the coherent and the in¬ 
coherent populations are in chaotic states (called as chaos-chaos chimeralike states) and, in another 
type, the incoherent population is in periodic state while the coherent population has irregular small 
oscillation. Interestingly, we also recorded a metastable state in a parameter regime of the Lienard 
system where the coherent and noncoherent states migrates from one to another population. To test 
the generality of the coupling configuration, we present another example of bistable system, the van 
der Pol-Duffing system where the chimeralike states are observed, however, the coherent population 
is periodic or quasiperiodic and the incoherent population is of chaotic in nature. Furthermore, we 
apply the coupling to a network of chaotic Rossler system where we find the chaos-chaos chimeralike 
states. 

PACS numbers: 05.45.Xt, 05.45.Gg 


I. INTRODUCTION 

Chimera states emerge [IHiQ] as sequentially organized 
subpopulations of coherent and incoherent dynamical 
units in a network of oscillators under nonlocal coupling. 
From first observation of this unexpected phenomenon 
in a network of phase oscillators BE] in the weak cou¬ 
pling regime, till date it has been reported to exist in 
limit cycle systems m m and chaotic systems m in 
the stronger coupling limit too. There in addition to 
phase incoherence, amplitude incoherence of a subpop¬ 
ulation has been found in the chimera states. Evidence 
of chimers states, by this time, has been found in chemi¬ 
cal El, opto-electronic [12] and electronic circuit experi¬ 
ments muni, and lately, in an experiment with network 
of metronomes El- 

Three different categories of chimera states have so far 
been identified m in networks of limit cycle or chaotic 
oscillators under nonlocal coupling. The basic chimera 
structure is composed of an incoherent subpopulation in 
a chaotic state while the coherent subpopulation could 
be periodic m 0 H 0 uni or remain close to a steady 
state m- In another type of chimera states [H [7] , the 
incoherent population remains in a state of spatial chaos 
El while the coherent population may be in a steady 
state or periodic state. A third kind of chimera state is 
classified as to coexisting structure of spatial chaos and 
spatio-temporal chaos [16] in the incoherent population. 
At least one bistable system was found m where all 
three types of chimera states exist in different parameter 
regimes, however, we emphasize that the network was 
under nonlocal coupling. 

Chimera states are intriguing since it emerges in an en¬ 
semble of identical oscillators under symmetric coupling 


although nonlocal. It is more nontrivial in an ensem¬ 
ble of identical oscillators under all-to-all global coupling 
since no spatial sequence or identity of the oscillators 
exists. However, a population of globally coupled oscil¬ 
lators was reported to split miiisiEa into synchronized 
and desynchronized subpopulations which has also been 
called as chimera states. We preferably call it chimeralike 
states as suggested by others [19] since there is no spa¬ 
tial pattern yet reminiscent of the chimera states under 
the nonlocal coupling. Such chimeralike states were no¬ 
ticed in the past EOlllIlllS] in globally coupled network, 
although not defined explicitly until stated categorically 
by Sen et al [22]. Almost at the same time, it has also 
been reported in limit cycle systems for a nonlinear global 
coupling m, globally coupled phase oscillators with de¬ 
layed feedback [19]. The mechanisms of the emergence 
of chimeralike states differ for different coupling configu¬ 
rations in limit cycle systems; it is either amplitude me¬ 
diated nil ns or amplitude modulated chimera [18] . In 
the chimeralike states too, the phase and/or amplitude 
of the coherent population are randomly distributed in 
the incoherent population and the coherent population is 
in periodic state. 

Nonisochronicity [22ll25] plays a crucial role in the 
chimeralike states of globally coupled network such as 
the case of Complex Ginzburg-Landau system [22] and 
the van der Pol system [26] . Otherwise a nonlinear global 
coupling m can also break the symmetry of a popula¬ 
tion into synchronous and nonsynchronous subpopula¬ 
tions. The presence of delay feedback as shown in a net¬ 
work of phase oscillators under global coupling may also 
create m such bistability of synchronous and nonsyn¬ 
chronous states in a population. Alternatively, a com¬ 
bination of attractive and repulsive coupling was also 
shown [19] to break the symmetry of globally coupled 
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bistable oscillators to create ctiimeralike states, however, 
they were forced into separate state variables of each dy¬ 
namical unit of the network from an external dynamical 
source. 

In contrast, in this paper, we use the attractive and 
repulsive global interactions in networks of oscillators, 
limit cycle and chaotic systems, to observe chimera- 
like states but do not apply the coupling interaction 
from an external dynamics m- We assume that the 
coupling interactions originate within the system. We 
use the attractive coupling as a mean-field self-feedback 
while the repulsive interaction is introduced either as 
a mean-field self-feedback or a cross-feedback. A type 
of self-feedback and cross-feedback coupling was used 
earlier [7] in a network under nonlocal coupling, where 
both attractive and repulsive interaction were present, 
to show chimera/multichimera states. We use a sim¬ 
ilar self-feedback as well as cross-feedback but purely 
global coupling. We first apply the coupling to a net¬ 
work of bistable Lienard system EH and find two types 
of chimeralike states in two different regions of the pa¬ 
rameter space of the network where the dynamics are 
qualitatively different in the synchronous and nonsyn- 
chronous populations. In one type, both the synchronous 
and asynchronous populations are chaotic, which is dif¬ 
ferent from the typically observed miiiniEa dynamics 
in the chimeralike states in limit cycle systems. In the 
second type of chimeralike states, the noncoherent pop¬ 
ulation is periodic but with no phase coherence which 
is qualitatively similar to the chimera states for nonlo¬ 
cal coupling reported earlier [7]. The coherent popula¬ 
tion shows a small oscillation close to the steady state 
but of irregular nature. Furthermore, we report various 
clustered states (1-, 2-, 3-, 4-cluster) and a special kind 
of dynamical behavior, namely, a metastable state in the 
network of the Lienard systems. In this metastable state, 
both the coherent and incoherent states migrate from one 
subpopulation to another in time, however, we find a dis¬ 
tinct network size effect in its transient behavior which we 
elaborate later. To further exemplify the role of the pro¬ 
posed coupling, we apply it to another network of bistable 
van der Pol-Duffing system and confirm the presence of 
chimeralike states where the noncoherent population is 
typically in chaotic state while the coherent subpopu¬ 
lation is periodic or quasiperiodic. Next for a network 
of chaotic Rossler systems, we simplify the coupling by 
separating the attractive self-feedback and the repulsive 
cross-feedback coupling from a single variable and apply 
both as self-feedback to two different variables and find 
clear evidence of the chimeralike states. The dynamics is 
typically chaotic in both the coherent and noncoherent 
subpopulations which we call as chaos-chaos chimeralike 
states. We elaborated the coupling structure in the next 
section, and we located the parameter regions of two dif¬ 
ferent types of chimerlike states, the clustered states and 
and the metastable state in the network of Lienard sys¬ 
tem in section III. The chimeralike states in the van der 
Pol-dufhng system and the Rossler system are described 


in sections IV and V respectively. Results are summa¬ 
rized in section VI. 


II. NETORK COUPLING CONFIGURTION 

The dynamics of the z-th node of the network is ex¬ 
pressed by, Xi = F(Xi,/i) -1- KAB where z = 1,A"; fi 
is the set of system parameters, K is the strength of cou¬ 
pling. All the dynamical nodes in the network are identi¬ 
cal, F : R^ ^ R^ (considering 2D systems here, extend¬ 
able easily to higher dimension); Xi = [xi^yi]"^ where z = 
F(Xi,/u) = A 

is a 2 X 2 matrix with real values and 5 is a 2 x 1 matrix 
defining two types of mean-field diffusions, 

V« 2 i a 22 ) \y-yij 

N N 

where ^ Xj and y= ^ ^ Vj- We discuss about 

3=1 i=i 

different options of the proposed coupling. Case 1: an = 
1 and ai 2 = a 2 i = a 22 = 0 describes the conventional 
global coupling, a type of self-feedback acting on a scalar 
variable x. Case II: an = 0, ai 2 = 0, a 2 i = e and 
a 22 = 1- The global coupling now consists of two com¬ 
ponents, one self-feedback involving the x-variable and 
another cross-feedback involving the ^-variable; they are 
both added to the dynamics of the ^-variable of each dy¬ 
namical unit of the network. Varying e from -\-ve to —ve 
value, the cross-feedback coupling changes from attrac¬ 
tive to repulsive nature. A combined effect of K and e on 
the collective and macroscopic behavior of the whole net¬ 
work is investigated. Case III: an = 1 and a 22 = e, other 
two elements are zero. The global coupling is established 
by applying two mean-field self-feedback interactions to 
the dynamics of two different state variables of each unit. 
Case IV: all the elements in matrix A are nonzero when 
it is of complex type. We focus here on the Cases II-III 
and explore chimeralike states in different example sys¬ 
tems. Note that we adopted a similar global coupling 
earlier [26] to observe the chimeralike states in the van 
der Pol system and the chaotic Rossler oscillator. We 
simplify the coupling here and show emergence of the 
chimeralike states in the bistable Lienard system, and a 
bistable van der Pol-Dufhng system and the Rossler sys¬ 
tem. Especially, we make elaborate studies on the net¬ 
work of the Lienard system and, locate regions of cluster¬ 
ing, two different kinds of chimeralike states and a type 
of metastable states in parameter space. 


III. CHIMERALIKE STATES: LIENARD SYSTEM 

We start our numerical study with an example of a 
Lienard system m and form a network of N number of 
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FIG. 1: Phase diagram in K — e plane. Different dynamical states are denoted as 1-cluster (dark red), 2-cluster (pink), 
3-cluster (red), 4-cluster (blue), chimeralike states (yellow). Various dynamics are shown at different parameter 
regimes: attractors of (a) Chimera-I {K = 1.6, e = —1.8), (b) 3-cluster {K = 1.7, e = —0.47), (c) metastable 
state {K = 0.91, e = —1.95), (d) 2-cluster {K = 1.7, e = —0.09), (e) 2-cluster {K = 0.7, e = —1.8), (f) 4-cluster 
(^K = 0.7, e = —0.59), (g) Chimera-II {K = 0.3, e = —1.8), time series of all the oscillators in (h) 1-cluster state 

(V = 0.2,6 = -0.2). 


identical units using the Case II coupling. 


A = 



fx-Xi\ 

\y-yij 


. The globally coupled network of the Lienard system is 

Xi = Vi, (1) 

yi = -axiyi - jSxf - ■yxi + K[{y - yi) + e(5 - Xi)](2) 

The Lienard system shows bistablity in isolation [23] : for 
a choice of parameters, a stable focus coexists with pe¬ 
riodic orbits. More categorically, for a choice of system 
parameters as given below, the system has a saddle sep- 
aratrix at (0,0) between a stable focus at (1,0) and a 
saddle focus at (—1,0). A homoclinic orbit (HO) exists 
at the saddle point (0, 0), which separates the state-space 
into two regions. For choice of initial conditions inside 


the HO, the system moves to the stable focus (1,0) after a 
transient; for choice of initial conditions outside the HO, 
multiple periodic orbits appear which have different fre¬ 
quencies, i.e, the system behaves like a non-isochronous 
system [23H25] . 

Next we draw a phase diagram to demarcate the pa¬ 
rameter regions of different macroscopic dynamics of 
the network, coherent or noncoherent states, two differ¬ 
ent chimeralike states using a statistical measure m, 
namely, a strength of incoherence (S). For this measure, 
the whole population is divided into M number of bins 
of equal length n = N/M and a local standard deviation 
cr(m) is then defined 


cr(m) = 




E 


{z)Y^ 




(3) 
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where m = 1,2, {z) = jj = Xi-Xi+i 

and < . >t defines a time average. Using this local stan¬ 
dard deviation, we measure the strength of incoherence 
(S) as 


5'=!- 


Z-^rn=l 




M 


Sm = Q{S - a{m)) 


(4) 


where 0(.) is the Heaviside step function, and ^ is a small 
predefined threshold. Chimera states and coherent states 
are distinguished by the S value, in general, for nonlo¬ 
cal and global coupling. The S = 0 identifies the in¬ 
coherent state while S = 1 defines the coherent state. 
The S lies in-between 0 < S' < 1 for chimera states. 
For global coupling as mentioned above, no spatial iden¬ 
tity or index exists for the dynamical nodes and hence 
the multichimera-like states, although appears, cannot 
be distinguished from the chimera states. What really 
matters in the chimeralike states in globally coupled net¬ 
work is the symmetry breaking of a population into co¬ 
herent and non-coherent subpopulations. Once the co¬ 
herent state is identified using the above statistical mea¬ 
sure, we use another algorithm to separate different clus¬ 
tered states: we record the instantaneous value of the Xi 
variable of all the oscillators after a long transient, and 
consider any two of them as belonging to a particular 
group wherever they are identical to each other within a 
small bound. Thus the oscillators having identical Xi val¬ 
ues forms a group; each group forms a separate cluster. 
Finally all such separate groups determine the number 
of clusters (l-,2-,3-,4- cluster). In a clustered state, the 
sum of the number of dynamical units in all the groups 
is the total number of units in the network. If this condi¬ 
tion fails we do not consider formation of the cluster. We 
further check the temporal evolution of all the oscillators 
for visual check. The metastable state is first identified, 
by the statistical measure, as the only noncoherent state 
in the parameter space. The dynamics in this regime is 
then looked into to recognize its unknown behavior as 
described below. 

Figure 1 shows the phase digram in — e parameter 
plane showing different dynamical regimes. We change 
K and e both in steps of 0.01 and use the fourth order 
Runge-Kutta algorithm to integrate the system with a 
time step size 0.01. The initial states for yi are chosen 
as yio = 2(1 - 4i) for 1 < i < ^ and yio = 2{^ - 3) for 
Y ^ ^ N with added small random fluctuations. 

All initial states for x-variables are set at zero. A black 
curve divides the K — e phase space into two regions. The 
network has trivial equilibrium points at (-1,0), 

(0,0) and (1,0). In the region below this black curve, the 
homogeneous steady state at (1,0) is stable while (—1,0) 
and (0, 0) are always unstable. Different clustered states 
coexist in this region as noted in the diagram. The emer¬ 
gent homogeneous steady state is consistent with earlier 
results [29] in a network under mixed attractive and re¬ 
pulsive global coupling. On the upper side of the black 
curve, the (1,0) fixed point becomes unstable besides 
the unstable fixed points (—1,0) and (0, 0) and there, all 


the states, clustered or chimera states, are robust to the 
choice of initial conditions. Different dynamical states 
are shown with their phase portraits in different regions 
of the parameter space : 1-cluster in yellow, 2-cluster in 
green, 3-cluster in red, 4-cluster states in dark red. In two 
different regions of 2-cluster states, the nature of the dy¬ 
namics are different although periodic. The chimeralike 
states are observed in the parameter space indicated by 
the blue regions. Above the black curve, the chimeralike 
states (Chimera-I in Fig. 1) are chaotic both in non¬ 
coherent and coherent populations which feature is un¬ 
common in limit cycle systems; most importantly, it is 
independent of the choice of initial conditions. Similarly, 
above the black curve, we find strips of 2-cluster (green), 
3-cluster (red)and 4-cluster states (dark red) indepen¬ 
dent of initial conditions. We find a noncoherent region 
(white) there which we identify later as metastable state. 
On the other hand, the chimeralike states (Fig. Ig) in 
the region below the black curve are limited to specific 
choices of initial conditions. There the dynamics of the 
noncoherent population is periodic when all the oscilla¬ 
tors are distributed in along the trajectory (open circles 
in the phase portrait). As a result the noncoherent pop¬ 
ulation has no amplitude and phase correlation while the 
coherent population is limited to small oscillations close 
to the steady state (1,0) but of irregular nature. Differ¬ 
ent clustered states (1-, 2-, 3-, and 4-cluster) are present 
below the black curve and coexist with the stable focus 
at (1,0). Note that no 1-cluster state exists above the 
black curve. 

To draw the stability line (black curve) that delin¬ 
eates the parameter space into two regions, we analyt¬ 
ically calculate the determinants and trace of the Jaco¬ 
bian matrix(J) at the trivial equilibrium points (—1,0), 
(0,0) and (1,0) in the K — e parameter plane. We find 
that the unstable foci ( — 1,0) and (0,0) remain unstable 
for the whole K — e phase space, only the stable focus 
(1,0) changes its stability when crossing the black line 
curve. The focus (1,0) becomes unstable when crosses 
the stability line to the upper region, the determinant 
J{det{J)) is 

det[J) = A 1 .A 2 ...A 7 V > 0 
and we calculate 

detiJ) = + (iV - l)|e](a - > 0 

so that, (-l)'^[a + {N - l)^e](a - > 0. 

But a + {N — l)^e = —1 and it implies 

(-ir(a-|e)^-i<0 (5) 

Ke > -1 (6) 

Hence the black curve in the K — e phase space is a 
rectangular hyperbola satisfying the equation Ke = —1 
(see APPENDIX for details). For Ke > —1, (1,0) is a 
stable fixed point below the black curve and for Ke < 
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— 1 it is unstable above. The black curve is numerically 
verified by calculating the eigenvalues of the Jacobian of 
the network at (1,0). It is interesting to note here that 
the black curve is independent of the network size and 
it indicates that the chimeralike states are independent 
of the network size. The Ke < —1 condition implies 
K and e are to be of opposite sign, i.e., a combination 
of attractive and repulsive coupling is necessarily to be 
chosen for the chimeralike states to emerge. 

Figure 2 presents a snapshot in polar coordinates 
and the spatio-temporal evolution of the x-variables 
for chimera-I state when K = 1.6 and e = —1.8. The 
synchronized population is clubbed into a square (red) 
shown in Fig. 2(a) while the circles (blue) represent 
the desynchronized population distributed in both 
amplitude and phase. Both the synchronized and 
desynchronized populations are in chaotic state but the 
attractors are separated in phase space as shown in their 
phase portrait in Fig. 1(a). Figure 2(b) represents the 
time evolution of the x^-variable of all the oscillators 
showing strips of coherent and noncoherent nodes which 
continue for a long run along the ^-axis. Multiple strips 
of noncoherent nodes are seen which should not be 
confused with multichimera states as explained above. 
Basically the whole population breaks into two coherent 
and noncoherent subpopulations with no positional 
identity. 


90 « 



(a) (b) 

FIG. 2: Chimera-I state in a network of Lienard system 
for K = 1.6 and e = —1.8. (a) Snapshot in polar co¬ 
ordinate of all the oscillators, synchronized population 
in square (red) and nonsynchronous population in cir¬ 
cles (blue), (b) temporal dynamics of all the nodes in the 
network. 

Figure 3 shows the chimera-II state for K = 0.3 and 
e = —1.8. Figure 3(a) shows a snapshot in polar coor¬ 
dinates where the synchronized population is again de¬ 
noted by a single square (red) and the desynchronized 
oscillators in circles (blue) as seen distributed in phase 
and amplitude. Figure 3(b) shows the temporal dynam¬ 
ics of all the nodes which clearly reveals the coherent and 
noncoherent population for a long time run. The attrac¬ 
tors of the synchronized and desynchronized regions are 


90 “ 




|2 


0 


-2 


Node index 


(a) (b) 

FIG. 3: Chimera-II state in the network of Lienard sys¬ 
tem for K = 0.3 and e = —1.8. (a) snapshot in polar 
coordinate of all the oscillators, synchronized population 
in square (red) and nonsynchronous population in cir¬ 
cle (blue), (b) temporal dynamics of all the nodes in the 
network. 


shown in Fig. 1(f). The desynchronized oscillators are 
all periodic, however, they do not have any phase coher¬ 
ence as seen in distributed circles denoting positions of 
the oscillators in 2D phase portrait. The dynamics of the 
synchronized population in a small dot inside the phase 
portrait is enlarged in the inset that shows small irregu¬ 
lar oscillation although remains very close to the steady 
state at (1,0). 

Most interestingly, an unknown kind of dynamics 
which we denote as metastable states appear near the 
boundary of the noncoherent state (white region) and 
the chimera-I state hi K — e phase space (Fig. 1). Both 
the coherent and noncoherent populations are chaotic 
and their attractors are overlapping each other in state 
space (Fig. Ic). In the perspective of cortical dynamics 
m, a relative coordination in neuron population has a 
temporal behavior; stronger coordination at one time 
and becomes weaker at another time and this relative 
coordination may switch within the population. Figure 
4 shows the time evolution of x-variable of all the 
oscillators in the metastable states for K = 0.91 and 
e = —1.95 for N=10 (left panels) and N=100 (right 
panels) respectively. The coherent population migrates 
randomly between the population in time as shown in 
the spatio-temporal plot in the left lower panel, which 
is defined as a metastable state. Goherence in the sense 
of coordination between a group of oscillators shows a 
temporal change. However, after a long transient, it 
moves to a cluster state as shown in the left upper panel 
when the number of oscillators is considered as N=10. 
For N=100, the metastable states clearly continues for 
a long time as seen in the right panels. We find that 
this state is a long lived transient and the transient time 
increases rapidly with the oscillator number. We plan to 
explore this metastable states, in further details, in the 
future. 
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(a) (b) 

FIG. 4: Temporal dynamics of Xi variable for all the 
oscillators for K = 0.91, e = —1.95. Clustering is seen in 
the left upper panel {N = 10) after a long time. Right 
panels (A^ = 100) show no clustering. Upper and lower 
panels are splitted in time to show the dynamics for long 
time in both the examples. 


IV. CHIMERALIKE STATES: VAN-DER 
POL-DUFFING SYSTEM. 

We construct a network of a bistable van der Pol- 
Duffing system m using the attractive self-feedback and 
repulsive cross-feedback global coupling {Case II). 

Xi = Vi, (7) 

yi = a{l — xl)yi — + Fsinwt 

+ K[e{x - Xi) + {y - yt)]. (8) 

and parameter values are a = 0.2, F = 1 and uj = 0.94 
when the isolated system is bistable having one peri¬ 
odic and one chaotic attractor. The initial states for 
Hi are chosen as i/io = 3(1 — |^) for i = 1 to y and 
= 3(1^ — 3) for i 1 V and the initial 

states for Xi are chosen as = 2(1 — |^) for i = 1 
to Y = 2(|^ — 3) for i ^ + 1 fo ^ with an 

added small random fluctuation. Figure 5 shows a snap¬ 
shot in polar coordinates (upper left) where the square 
(red) represents the synchronized population and the dis¬ 
tributed circles (blue) correspond to the desynchronized 
oscillators. In the chimeralike states, the synchronized 
population is in periodic state (could be quasiperiodic 
for a different choice of parameters) and the desynchro¬ 
nized oscillators are in chaotic states. Right upper panel 
shows the time evolution of x-variables of all the oscil¬ 
lators which confirms the coexistence of coherent and 
noncoherent subpopulations and the number of oscilla¬ 
tors in each population remain unchanged for a long run. 
Again we ignore the concept of multichimera for reasons 
explained above although the typical signature of multi¬ 
chimera as found for nonlocal coupling exists here. Lower 
panels again show two clustered state for K = OA and 


e = —4 when the left panel clearly identifies two ciustered 
popuiations in a circie (bine) and a square (red) in poiar 
coordinate and time evoiution piots of aii the osciiiators. 
A smaii change in the A-vaiue originates the chimeraiike 
states. 



(a) 


(b) 
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FIG. 5: Ghimera state for K = 0.2 and e = —4 and two 
ciuster state for K = OA and e = —4 in van der-Poi- 
Dufhng system. Poiar piots at ieft, upper panei shows 
random distribution of phase and ampiitude of the non- 
herent osciiiators in circies (bine) and singie coherent 
ciuster in square (red). Right paneis piot spatio-temporai 
dynamics. Upper panei shows chimeraiike state when aii 
the nodes spiits into coexisting coherent and noncoherent 
subpopuiations. Lower right panei shows two ciustered 
popuiations. 


It is worth mentioning that both the attractive and 
the repuisive mean-fieid coupiing may be appiied as seif- 
feedback, in the sense, that they are added separateiy to 
their corresponding state variabies; the cross-giobai feed¬ 
back is not a necessary condition to observe chimeraiike 
states in the exampie systems. We do not present the 
resuits here for the systems, discussed above, however, 
eiaborate this using the chaotic Rossier modei. 


V. CHIMERALIKE STATES: ROSSLER SYSTEM 

We simpiify the coupiing [26] by separating the attrac¬ 
tive and the repuisive mean-fieid interactions and appiy 
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(a) (b) 

FIG. 6: Chimera state for Rossler oscillators for Ki = 
0.071 and K 2 = 0.151. Snapshot of Xi{t) at left and 
phase Oi at right for all the oscillators. 


them as self-feedback to two different dynamical equa¬ 
tions {Case III). We choose a network of Rossler oscilla¬ 
tors 


Xi = -Vi - Zi - Ki{x - Xi) (9) 

yi = Xi + ayi + K2{y - Vi) ( 10 ) 

Zi = hxi + Zi{xi - c) ( 11 ) 

and the parameter values as a = 0.36, b = 0.4 and 


c = 4.5 in the chaotic regime and here the isolated 
system is not bistable. Figure 6 shows snapshots of 
x^-variable and the phase <pi for N = 100 oscillators. 
Left panel is a snapshot of amplitude and the right shows 
snapshot of phase of all the oscillators (nodes) showing 
signature of chimeralike states. The distribution of 
phase and amplitude along the nodes confirms the state 
of incoherence in one subpopulation while the other 
counterpart remain synchronized. We emphasize once 
again that it should not be confused with multichimera 
states. The whole population splits into coherent and 
noncoherent subpopulations with no specific spatial 
structure. 


VI. CONCLUSION 

We observed chimeralike states in networks of iden¬ 
tical nonlinear oscillators using a global coupling con¬ 
sisting of both attractive and repulsive mean-field feed¬ 
back. Historically, the chimera states were observed [T], 
as a strange phenomenon, in network of identical os¬ 
cillators under range limited interaction, the so called 
nonlocal coupling. The chimera states appeared highly 
nontrivial in globally coupled identical oscillators, how¬ 
ever, the chimeralike states were reported in globally cou¬ 
pled network of oscillators. A homogeneous network of 
globally coupled oscillators splits into coexisting coher¬ 
ent and noncoherent subpopulations in a selected pa¬ 
rameter space and for strong coupling. The coupling 
scheme could be established as a simple global mean- 
field interaction but the presence of nonisochronicity was 


found crucial for chimera states to observe [22l [26] in 
the limit cycle systems such as the Complex Ginzburg- 
Landau system, the van der Pol system. Alternatively, a 
nonlinear coupling [18] or a delay feedback m was used 
for chimeralike states to emerge in globally coupled net¬ 
work. A bistabilty in the dynamical units was found to 
augment the emergence of chimera states in limit cycle 
systems. However, for chaotic systems, neither the non¬ 
isochronicity nor bistability criterion is a necessary for 
the origin of chimeralike states as first shown in globally 
coupled chaotic map [21] . We showed here that a combi¬ 
nation of attractive and repulsive mean-field coupling can 
produce chimeralike states in a broad parameter space of 
a Lienard system. Interestingly, we identify two different 
types of chimeralike states. In one type, we found both 
the subpopulations in chaotic state. In another type of 
chimeralike states, the coherent population remains close 
to a steady state but with a small irregular oscillation. 
On the other hand, the noncoherent population remains 
in periodic or quasiperiodic state. Thus, in the globally 
coupled network under two competitive coupling compo¬ 
nents, a rich variety of chimeralike states with diverse dy¬ 
namical features were found. We established the role of 
two competing coupling components in creating chimera¬ 
like states using numerical examples of two bistable limit 
cycle systems, namely, a Lienard system and the van der 
Pol-Duffing system, and the chaotic Rossler system. We 
noted that such attractive and repulsive coupling were 
used m in globally coupled chaotic systems recently to 
observe chimeralike states but they were forced into the 
network from an additional external dynamics. In con¬ 
trast, we explored the chimeralike states in both limit 
cycle and chaotic systems where the coupling interac¬ 
tions are all internally generated and applied either as 
a self-feedback and/or cross feedback. To make a clear 
distinction of the chimera states in globally coupled net¬ 
work from the nonlocally coupled network, we prefer the 
terminology such as the chimeralike states throughout 
the text. This is due to the fact that chimeralike states 
refers to a phenomenon of symmetry breaking of a homo¬ 
geneous population into subpopulations of coherent and 
noncoherent oscillators; no spatial identity exists. On 
the other hand, in the traditional chimera states, a clear 
spatial identity of oscillators exists due to the nonlocal 
nature of the coupling. The whole population splits into 
coherent and noncoherent oscillators but organized in a 
spatial order. 

Besides chimeralike states, we observed different clus¬ 
tered states (1-, 2-, 3-, 4-cluster) and, most importantly, 
a kind of metastable state in a parameter region near the 
transition from a cluster state to the chimeralike states 
via a noncoherent state. In the metastable states, the co¬ 
herent/noncoherent population migrates in time to dif¬ 
ferent population of oscillators. At one instant, some 
oscillators were in coherent state which loses coherence 
with time but another population establishes coherence 
by that time. This migrating coherent or noncoherent 
state led to a permanent coherent state after a long tran- 
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sient for network of smaller size {N = 10). But this 
transient behavior showed a network size effect. For a 
reasonably larger network size {N = 100), we could not 
predict the transient time of the coherence state even for 
large simulations within our limited computational facil¬ 
ity. We plan to explore the metastable state, in further 
detail, in the future. 

S.K.D. and P.K.R. acknowledge individual supports by 
the CSIR Emeritus Scientist Scheme, India. A. Mishra 
is supported by the University Grant Commission, India. 


the rows where only one element is 1 and all other ele¬ 
ments are zero.) and the reduced determinant is a A" x A" 
determinant and it is expressed by 



K 

K 

K 

a 


••• 

••• 






a 









a 



det{J) = (—1)^ 


APPENDIX 

For coupled Lienard system the Jacobian matrix 
at(l, 0) can be written as 

J = 1^U + -V) (12) 




In the first step we add the remaining rows to the first 
row and then pull out constant out of the determinant. 

det{J) = 


0 0 
K 

N 


with a = —3yd — 7 + — Ke and b = —a ^ — K. 

1 is the unit matrix of order A x A and I is the identity 
matrix of order A x A. 

So the Jacobian matrix(J) at (1,0) becomes a 2A x 2A 
square matrix and it is denoted as 


1 1 1 . 1 


a 


a 


-l)^[a+{N-l)-e] 


1 

0 

0 

0 

0 ••• 

... 0 

O' 


h 


K 


K 

... ^ f 

K 


U 


N 


N 


N 


0 

0 

1 

0 

0 ••• 

... 0 

0 

In the next 

Ri ^ Ri — ( 

K 

N 

a 

b 

N ^ 

K 

N 

... Ae 

K 

N 

i = 2,3,...A. 

0 

0 

0 

0 

1 ••• 

... 0 

0 

det{J) = 


\K K . K 

\ a I 

In the next step we perform the row operation 
Ri Ri — where Ri denotes i-th row with 


0 (a - f e) 


K K K K K K 
-N^ N N N 


(-l)^[a+(7V-l)-e] 


For (1,0) to be stable fixed point, two conditions must 
be satisfied 


det[J) = A 1 .A 2 ...A 7 V > 0 
Tr{J) = Ai -F A 2 + ... + Atv < 0 


(a-fe) 


with Ai, A 2 ,..., Aat are the eigenvalues of J. 

First we use Laplace expansion of the determinant of 
J{det{J)) with respect to the odd-numbered rows (i.e. 


det{J) = (-l)^[a + (AT - l)-e]{a - (15) 
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or 

(-l)^[a+(AT-l)^e](a-> 0 (16) 

Now putting value of a, we get a + {N — l)^e = —1. 
This implies that 

{a-^e)^-^<0 (17) 

Ke > -1 (18) 

Hence the solid black curve in the K — e phase space is a 
rectangular hyperbola having the equation Ke = —1. For 


Ke > —1, (1,0) is a stable fixed point and for Ke < —1, 
(1, 0) is unstable. 

Again from equation(14) 


Tr{J) < 0 

(19) 

bN <0 

(20) 

Na 

(21) 



Equations (18) and (21) must be satisfied for (1,0) to be 
a stable fixed point in the network. Similarly we check 
the stability of ( — 1,0) and find that it is unstable in the 
whole K — e phase space. 
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